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FOR  SYSTEM  RELIABILITY 

by 

M.  Vernon  Johns,  Jr. 

ABSTRACT 

The  basic  problem  of  determining  objective  ( frequent istic) 
confidence  bounds  for  the  reliability  of  a  series  system  based  on 
failure  data  from  tests  of  the  independent  components  is  addressed. 
The  notion  of  confidence  bounds  based  on  orderings  imposed  on  the 
sample  space  is  exploited,  and  certain  optimality  considerations  are 
incorporated.  Advantage  is  taken  of  the  simplifications  resulting 
from  the  use  of  the  Poisson  approximation  for  data  from  highly 
reliable  components.  Tables  of  exact  confidence  bounds  are  produced 
for  the  case  of  two-component  systems.  These  bounds  are  computed 
using  sample  orderings  generated  sequentially  by  a  two-stage,  pro¬ 
spective  optimization  procedure.  A  generalization  of  the  Lindstrom- 
Madden  technique  is  proposed  for  using  the  tables  to  find  confidence 
bounds  for  systems  consisting  of  more  than  two  components  with 
differing  sample  sizes. 

Key  Words:  Reliability,  series-system,  confidence  bounds,  Poisson 
approximation,  Lindstrom-Madden,  sample  orderings. 
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FEASIBLE  OBJECTIVE  CONFIDENCE  BOUNDS 
FOR  SYSTEM  RELIABILITY 

By 

M.  Vernon  Johns,  Jr. 

1.  Introduction 

Certainly  one  of  the  most  basic  statistical  problems  in  the 
assessment  of  system  reliability  is  that  of  determining  a  confidence 
bound  on  the  reliability  of  a  series  system  based  on  component  fail¬ 
ure  data.  The  continuing  appearance  of  papers  concerned  with  this 
subject  (e.g.,  Harris  and  Soms  1981,  Butcher  et  al.  1978,  and 
Winterbottom  1980)  testifies  to  the  elusiveness  of  solutions  which 
are  both  practically  feasible  and  acceptably  precise.  The  present 
paper  deals  with  the  case  of  systems  characterized  by  high  intrinsic 
reliability  (>90%)  where  the  use  of  the  Poisson  approximation  for  the 
binomial  distributions  of  component  failure  data  does  not  introduce 
appreciable  error.  The  emphasis  is  on  objective,  frequentistic  con¬ 
fidence  bounds  which  avoid  the  uncertainties  of  interpretation  asso¬ 
ciated  with  posterior  bounds  obtained  by  Bayesian  methods.  The  use 
of  the  Poisson  distribution  provides  several  advantages.  It  easily 
leads  to  valid  results  for  cases  involving  zero  observed  failures  for 
some  components  where  maximum  likelihood  methods  and  other  proposed 
approximations  tend  to  break  down  (see,  e.g.,  Madansky  1965  and  Mann 
et  al.  1974).  It  also  permits  the  pooling  of  failure  data  for  dif¬ 
ferent  components  having  the  same  test  sample  sizes.  This  potential 
reduction  in  the  effective  number  of  system  components  enhances  the 
utility  of  tabulated  bounds  such  as  those  presented  here. 
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Because  of  the  structure  of  the  problem,  universally  optimal 
confidence  bounds  for  system  reliability  (i.e.,  "uniformly  most 
accurate"  bounds  in  the  sense  of  Lehmann  1959)  do  not  generally  exist 
On  the  other  hand,  the  ideas  of  Buehler  (1957)  may  be  exploited  to 
produce  a  variety  of  valid  confidence  bounds  based  on  total  orderings 
of  the  sample  points.  Such  bounds  are  exact  in  the  sense  that  the 
desired  coverage  probability  is  guaranteed.  The  construction  of  good 
confidence  bounds  is  thus  reduced  to  the  selection  of  suitable  order¬ 
ings  imposed  on  the  sample  space.  This  is  the  approach  adopted  in 
the  present  study.  Previous  applications  of  these  methods  to  relia¬ 
bility  may  be  found,  for  example,  in  Harris  and  Soms  (1980),  Johnson 
(1969) ,  and  Lipow  and  Riley  (1959) . 

Once  the  sample  ordering  approach  has  been  chosen,  there 
remains  the  problem  of  determining  orderings  which  lead  to  confidence 
bounds  which  are  "good"  according  to  some  criterion  measuring  the 
size  (length)  of  the  confidence  region.  An  early  proposal  of  the 
present  author  (Johns  1975)  was  to  generate  the  ordering  by  means  of 
a  simple  function  of  the  observations  which  was  asymptotically  equiv¬ 
alent  to  the  maximum  likelihood  confidence  bound.  This  method  guar¬ 
antees  asymptotic  bptimality  when  the  numbers  of  component  failures 
observed  under  testing  is  large.  It  was  found,  however,  that  the 
resulting  bounds  could  be  noticeably  improved  for  small  to  moderate 
numbers  of  observed  failures  by  more  sophisticated  methods.  Another 
procedure  investigated  (Johns  1977)  was  a  sequential  method  for  gen¬ 
erating  the  sample  ordering  starting  at  the  origin  (zero  component 
failures  observed)  and  selecting  at  each  stage  as  the  next  sample 
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point  the  "adjacent"  point  producing  the  largest  value  for  the  lower 
hound  on  reliability.  This  method,  while  intuitively  appealing,  does 
not  generally  produce  a  "best"  ordering  as  has  been  suggested  by  some 
investigators.  In  particular,  it  is  improved  upon  by  the  method 
adopted  in  the  present  paper. 

The  non-existence  of  a  unique  optimal  ordering  leaves  open 
the  possibility  of  obtaining  a  confidence  bound  which  is  at  least 
admissible  by  choosing  the  sample  ordering  to  minimize  the  expected 
length  of  the  confidence  interval  under  some  suitable  prior 
distribution.  Such  a  semi-Bayesian  approach  does  not  in  any  way 
impair  the  frequentistic  interpretation  of  the  confidence  bounds 
obtained  from  the  minimizing  ordering.  The  implementation  of  such  a 
minimization,  while  theoretically  perfectly  possible,  turns  out  to  be 
totally  unfeasible  computationally  except  for  the  earliest  part  of 
the  ordering  generated.  Nevertheless,  for  a  class  of  priors  chosen 
to  emulate  certain  properties  of  maximum  likelihood,  fragmentary 
orderings  computed  by  this  method  provide  a  considerable  justifica¬ 
tion  for  the  two-stage  "look-ahead"  sequential  method  actually  used  to 
generate  the  tables  which  are  a  principal  concern  of  this  paper. 

Suppose  that  the  series  system  under  consideration  consists 
of  k  independent  components  and  that  the  respective  probabilities  of 
component  failure  are  p±  =  1  -  q  i  =  l,2,...,k.  The  system  relia¬ 
bility  R  is  given  by 

k  k 

r=  n  q  -  n  (i-p)  .  (l.D 

i=l  i=l  1 
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If  the  observed  numbers  of  failures  for  the  k  components  are 

based  on  independent  tests  with  corresponding  sample 
sizes  n^,n2, . . .,nfc,  we  let  =  n^p^  i  =  l,2,...,k,  so  that 


R  =  II  (1-  X./n.)  =  1  -  E  X./n, 
i=l  1  i=l  1  1 


(1.2) 


The  approximation  on  the  right  will  be  best  when  the  p^'s  (=  X^/n.) 
are  all  small  which  is  just  the  case  where  R  is  close  to  one  and  the 
Poisson  approximation  for  the  distributions  of  the  X^’s  is  valid. 

It  will  be  convenient  to  express  the  problem  in  a  canonical 
form  by  introducing  some  further  notation.  Let  c  =  E^_^  1/n.  and 
=  1/cn^  i  =  l,2,...,k.  Then  letting  X  =  (X^.X^  . . .  ,Xk)  and 


k 

0(X)  =  E  a  X  ,  (1.3) 

i=l  1  1 


we  have 


k 

R  =  1  -  c  E  a  X.  =  1  -  c0(X)  .  (1.4) 

i=l  1 

We  shall  assume  henceforth  that  the  components  are  indexed  so  that 
_>  >  . . .  _>  n^  which  implies  that  a^  £  ±  . . .  <_  a^.  The  pro¬ 

blem  of  finding  a  lower  confidence  bound  for  R  is  thus  reduced  to 
that  of  finding  an  upper  confidence  bound  for  0 (X)  .  The  fact  that 
0 (X)  is  a  convex  combination  of  the  X^s  facilitates  tabulation  of 
the  bounds  by  reducing  the  number  of  classification  variables  by  one. 
In  principle,  confidence  bounds  for  R  could  be  constructed  directly 
without  introducing  the  approximation  (1.2).  Such  an  approach  would. 
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however,  eliminate  the  possibility  of  constructing  useful  tables  of 
bounds,  since  separate  entries  would  be  required  for  every  configura¬ 
tion  of  (n^n^,  . .  .,1^)  . 

The  vector  of  observations  is  X  =  (X^X^  . . .  ,Xk)  so  that  the 
sample  space  on  which  a  total  ordering  must  be  imposed  consists  of 
all  vectors  x  =  (x^x^ . . . ,xk)  where  the  x^’s  are  non-negative 
integers.  It  seems  reasonable  to  confine  attention  to  orderings 
which  are  consistent  with  the  natural  partial  ordering  induced  by 
dominence  (see  Section  2),  and  we  shall  do  so.  In  Section  2  it  is 
shown  that  the  best  upper  confidence  bound  for  0(A)  with  confidence 
coefficient  1  -  a  which  is  monotone  in  a  prescribed  total  ordering 
(designated  by  the  relation  ^  )  is  given  by 

t(x)  =  sup  0(A)  ,  (1.5) 

AeS*(x) 

or 

where  S*(x)  =  (A  :  P,  {X  x}  =  a)  .  From  (1.4)  the  lower  confidence 

CX  A 

bound  for  system  reliability  is  then  given  by 

r(x)  =  1  -  ct(x)  .  (1.6) 
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Because  of  the  additive  property  of  independent  Poisson 
observations  the  effective  number  of  components  in  the  system 
may  be  reduced  if  some  sample  sizes  are  equal  (or  nearly  equal)  as 
follows:  Suppose  n^  =  n^  =  . ..  =  n  «  n  for  some  r,  2  <  r  <  k.  Then 
a.  =  a_  =  ...  =  a  =  a*  (say)  and  letting  A*  =  E1;  ,  A.  we  have 

X  IT  1 


0(A)  =  a*A*  + 


k 

E 

i=r+l 


aiXi 


(1.7) 


Now  X  =  E^_^  has  a  Poisson  distribution  with  parameter  A  so  that 
the  right  hand  side  of  (1.7)  has  the  correct  form  for  the  case  of 
dimension  k*  =  k-r+1  except  that  the  coefficients  must  each  be 
divided  by  c^  =  a  +  avJ/|  +  ...  +  at  to  preserve  convexity .  The  con- 


r+1 


fidence  bound  for  the  k*  dimentional  case  based  on  X*,Xr+^, . . . ,X^  and 
the  normalized  coefficients  may  then  be  multipled  by  c^  to  obtain  the 
bound  for  the  original  9(A)  given  by  (1.7).  Further  reductions  may 
be  made  in  the  same  way  if  several  groups  of  components  have  common 
sample  sizes.  If  n^,^, . . . ,n^  are  only  approximately  equal,  their 
average  may  be  used  for  n  in  the  above  calculations  to  obtain  an 
approximate  bound. 

If  all  components  are  subjected  to  the  same  number  of  trials, 
we  may  take  k  equal  to  one  and  the  problem  reduces  to  the  familiar 
one  of  finding  an  upper  confidence  bound  for  a  single  Poisson 
parameter . 

The  component  failure  data  may  be  developed  through  independ¬ 
ent  testing  of  the  components,  or  through  testing  of  the  complete 
system  with  the  assignment  of  failures  to  the  appropriate  components. 
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Even  in  the  latter  case  component  sample  sizes  may  differ  if  com¬ 
ponents  are  redesigned  during  the  course  of  testing  so  that  the 
trials  and  failures  observed  prior  to  redesign  are  not  relevant  to 
the  reliability  of  the  final  version  of  the  system. 

It  is  worth  noting  that  the  same  formal  confidence  bound 
problem  arises  from  a  somewhat  different  and  more  restrictive  model 
involving  continuous  time  sampling  and  exponentially  distributed 
failure  times.  Thus,  if  we  assume  that  the  time  to  failure  of  the 
i-th  component  has  the  exponential  probability  density  function 

"V 

VU  e  ,  >  0,  t  >  0,  and  if  the  fixed  aggregate  testing  time  for 

the  i-th  component  is  given  by  T\  and  the  total  number  of  failures 
observed  is  X^,  i=l,2,...,k,  then  the  X^'s  which  are  sufficient  sta¬ 
tistics  have  independent  Poisson  distributions  with  corresponding 
parameters  =  y^x^.  Such  test  data  would  be  generated,  for 
example,  if  a  single  unit  of  each  component  were  placed  on  test  and 
repaired  whenever  a  failure  occurred  until  the  total  test  time  for 
each  component  reached  its  corresponding  limit,  x^.  The  system  reli¬ 
ability  for  a  single  duty  cycle  (consisting  of  operation  for  a  unit 
interval  of  time)  is  given  by 

k  -y .  k  1  k 

R  =  II  e  1  =  exp  -  II  y.  =  exp--  2  A  /x .  ,  (1.8) 

i-i  l  i-i  XJ  [  i-i  1  4J  . 
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and  the  confidence  bound  problem  is  essentially  the  same  as  the  one 
previously  introduced  except  that  no  approximations  are  required. 

The  general  theory  of  confidence  bounds  based  on  sample 
orderings  is  discussed  in  Section  2.  In  Section  3  the  case  of  sys¬ 
tems  having  two  components  is  considered  in  detail,  and  the  rational 
for,  and  use  of,  the  tables  for  this  case  are  explained.  Section  4 
includes  suggestions  for  constructing  approximate  confidence  bounds 
for  cases  of  systems  with  k  _>  3  by  finding  approximately  equivalent 
cases  with  k  =  2.  Use  of  the  maximum  likelihood  ratio  bounds  for 
cases  where  the  data  are  beyond  the  limits  of  the  available  tables  is 
also  discussed . 

2.  Bounds  and  Orderings:  Generalities 

The  idea  of  using  sample  orderings  to  generate  confidence 
bounds  was  first  introduced  by  Buehler  (1957)  who  discussed  the 
validity  of  the  proposed  method  in  the  context  of  a  specific  relia¬ 
bility  problem.  Bol'shev  and  Loginov  (1969)  discuss  the  construction 
of  confidence  bounds  monotone  in  the  sample  orderings  generated  by 
certain  functions  of  the  observations.  In  the  following,  which  is  a 
revision  and  extension  of  Johns  (1975) ,  we  develop  the  theory  with 
emphasis  on  the  sample  orderings  themselves  rather  than  possible 
generators  of  the  orderings. 
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To  develop  the  general  ideas  relating  exact  confidence  bounds 
to  sample  orderings  it  is  convenient  to  introduce  a  fairly  abstract 
statistical  model.  Let  the  sample  space  X  he  endowed  with  a  measura¬ 
ble  total  ordering  relation  "  ^  "  and  let  X  represent  the  random  out¬ 
come  of  the  experiment  where  the  space  of  possible  outcomes  is  I. 
Suppose  that  the  possible  distributions  of  X  are  determined  by  the 
family  of  probability  measures  P^,  indexed  by  A,  an  element  of  the 
parameter  space  A.  Our  objective  is  to  find  a  1  -  a  level  upper  con¬ 
fidence  bound  for  a  specified  real-valued  function  0(A)  defined  on  A 
where  the  range  T  of  0(A)  is  assumed  to  be  closed  and  bounded  below. 
The  quantity  ae  (0,1)  is  regarded  as  fixed  throughout.  We  make  the 
following  definitions  and  assumptions: 

Definition  Dl.  For  each  x  e  X  let 

S(i(x)  =  (A  :  P^(X^.  x}  >  a}  .  (2.1) 


Definition  D2.  For  each  xeX  let 


t(x)  = 


sup{0(A)  :  A e  Sa(x)}  ,  if  Sa(x)  is  non-empty  , 


[inf  T,  otherwise 


(2.2) 


Remark  1.  By  D2  if  0(A)  >  t(x),  then  necessarily  P^{x^.  x)  £  a. 


Assumption  Al.  For  every  subset  C  of  X  having  the  property  that  if 

x  e  C  and  y  ^  x  then  y  £  C,  there  exists  an  ordered  sequence 

X,  ^  x  ^  . . .  of  elements  of  C  such  that  C  =  lT  ,  (x  :  x*<  x  }  . 

n=x  '  n 

Assumption  A2.  For  each  x  £X,  if  0(A)  =  t(x),  then  P^{x<j.x}  £  a. 
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Remark  2.  By  Dl,  D2,  and  A2  if  0(X)  =  t(x),  then  X  t  S„(x),  i.e., 

vX 

the  supremum  in  D2  is  never  attained. 

We  now  establish  the  following  propositions. 

Proposition  Pi.  The  function  t(x)  is  monotone  in  the  ordering  on  X» 

Proof ;  If  x,yeX  and  x^  y,  then  S^(x)  a  Sa(y)  (Dl)  which  in  turn 
implies  t(x)  t(y)  (D2)  .  fj 

Proposition  P2 .  Under  assumptions  A1  and  A2  the  function  t(x)  is  an 
upper  confidence  bound  for  9(X)  at  level  1  -  a.  In  particular, 

PA(0(X)  <  t(X)}  >  1  -  a  for  all  Xe  A  .  (2.3) 

Proof :  For  arbitrary  XQ  eA,  let  0Q  =  9(AQ)  and  CQ  =  (x  :  t(x)<^0Q}. 

The  result  follows  immediately  for  all  Xq  for  which  is  empty. 

Assume  that  Cq  is  non-empty.  Then  by  PI  the  set  Cq  possesses  the 

property  required  in  A1  for  the  existence  of  a  sequence  {x^}  c:  Cq 

00  r  i 

such  that  x  ^  x  ...  for  all  n,  and  C_  =  U  -  {x  :  x^(  x  }.  This 

n  ^  n+1  ’  0  n=l  ”  n 

implies  that,  as  n  00 , 

P,  {X<x  }  *  P.  (Xe  Cn}  .  (2.4) 

Xo  n  X0  0 

But  by  Remark  1  and  A2,  for  all  n,  P^  (x^,  x^}  <  a.  Hence 
P^  (Xe  CQ}  a  and  the  desired  result  follows.  £] 
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Proposition  P3.  Under  assumptions  A1  and  A2,  if  t(x)  is  any  T-valued 
confidence  bound  such  that  P^{0(X)  <  t(X)}  >_  1  -  a  for  all  X e A,  then 

(i)  sup  t(y)  t(x)  for  all  x  el  ,  and 
y 

(ii)  if  t(x)  is  monotone  in  the  ordering  on  X,  then 
t(x)  t(x)  for  all  xeX  . 

Proof :  First  we  assume  that  t(x)  is  monotone  and  establish  (ii). 

Suppose  there  exists  an  x'  eX  such  that  t(x')  <  t(x')  •  Then  S  (x') 

a 

must  be  non-empty  and  by  Remark  2  following  A2,  the  sup  defining 

t(x')  is  not  attained.  Hence  there  exists  a  X'  eS  (x')  such  that 

a 

t(x')  <  0 (X  * )  <  t(x'),  and  P^tfX^.x’}  >  a.  Thus  by  the  monotonicty 
of  t(x) , 

Px,{t(X)  <  0(X*)>  >  Px,{t(X)  <  t(x’)}  -  PX,(X<  x'}  >  a  .  (2.5) 

v 

This  contradicts  the  hypothesis  that  Px(0(X)  <  t(X)}  >_  1  -  a  for  all 
X  and  establishes  (ii) .  To  show  (i)  we  let  t*(x)  =  sup  t(y)  and 
observe  that  t*(x)  is  monotone  in  the  ordering  on  X  and  t*(x)  t(x) 
for  all  x  eX-  Hence  if  t(x)  is  a  1  -  a  confidence  bound  for  0(X),  so 
is  t*(x)  and  applying  (ii)  to  t*(x)  yields  (i)  .  [] . 

In  order  to  specialize  these  results  in  the  direction  of 
applications  we  henceforth  assume  that  the  parameter  X  and  the 


11 


observation  X  are  both  of  dimension  k,  i.e.,  A  =  (A  ,A  ,  ...,A,)  and 

Z  K 

X  =  (X^,X2» . . • ,  X^)  where  the  A^'s  are  real  and  the  X^’s  are  random 
variables.  Without  essential  loss  of  generality  we  assume  that  A 
contains  the  positive  orthant.  Within  this  framework  we  make  the 
following  additional  assumptions: 


Assumption  A3.  The  function  0(A)  is  continuous  and  strictly 
increasing  in  each  of  the  A^'s. 

Assumption  A4.  For  any  xeX,  P-^Cx^x}  is  continuous  in  each  of  the 


Proposition  P4.  Assumptions  A3  and  A4  imply  that  Assumption  A2  is 
satisfied. 

Proof:  Suppose  that  for  some  x'el  there  exists  a  A'eA  such  that 

0 (A  )  =  t(x  )  and  P^t{x=$  x'}  >  a.  Then  by  A3  and  A4  we  can  find  a 

A  e A  with  A.  >  A.  for  all  i  and  A.  >  A.  for  some  i  such  that 
i  —  i  i  i  o 

o  o 

0(A")  >  t(x')  and  P^.i{x^x'}  >  a.  This  contradicts  D2  and  the 
result  follows.  [J 


Corollary  Cl.  If  (i)  the  observations  X^X^...^  are  independent 
Poisson  random  variables  with  parameters  A^,A2,...,Ak  respectively, 
and  (ii)  0(A)  =  a.jA^  +  a^2  +  •••  +  w^ere  the  a^'s  are  positive, 

then  t(x)  given  by  D2  is  an  upper  confidence  bound  for  0(A)  at  level 
1  -  a. 
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Proof:  A1  is  satisfied  for  any  ordering  since  X  is  discrete.  A3  is 
clearly  satisfied  for  0(A)  of  the  form  given  by  (ii) ,  and  A4  is 
satisfied  for  Poisson  random  variables.  The  desired  result  follows 
from  P4  and  P2.  j] 

The  actual  computation  of  the  bound  t(x)  given  by  D2  is 
greatly  facilitated  if  P^{x^  x)  is  monotone  in  the  components  of  X. 
The  following  proposition  gives  conditions  guaranteeing  this 
property: 

Assumption  A5.  The  components  of  X  are  independent  and  for 
i  =  1» 2, . . . ,k  the  distribution  of  each  X^  depends  only  on  the 
corresponding  A^^  and  is  stochastically  increasing  in  A.. 

Definition  D3.  We  denote  by  (  )  the  natural  partial  ordering  of  I 

generated  by  componentwise  dominance.  That  is,  x^y  if  and  only  if 

* 

y^  for  i=l,2,...,k,  with  x-^y  if  at  least  one  of  these 

inequalities  is  strict.  An  arbitrary  total  ordering  ( ^  )  on  X  is 
,  *  * 

consistent  with  (  )  if  x^y  implies  x  ^  y. 

Proposition  P5.  If  the  ordering  on  X  is  consistent  with  the  natural 
partial  ordering  and  Assumption  A5  is  satisfied,  then  for 
each  xeX,  P^fx-^x)  is  non- increasing  in  each  component  of  A. 

Proof :  For  y  e  X  we  introduce  the  representation  y  =  (y^,y^^)  where 
(2)  _  . 

y  -  (y2»y3»---»yk)-  For  real  z,  fixed  xeX,  and  all  yeX,  let 
Ix(z»y(2))  be  the  indicator  function  of  the  set  {y^  :  x)  . 
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Then 


P,{x4xl  -  E,  I  (X,,X(2))  .  (2.6) 

A  A  x  1 

For  any  yel,  if  z'  <  z",  then  (z',y^)  (z",y^),  by  the  con- 

(2)  (2) 

sistency  hypothesis,  and  I  (z',y  )  _>  I  (z",y  ).  Hence  letting 

X  X 

/n\ 

G  (z)  =  E.{l  (X  ,X  ;) (Xn  =  z}  we  see  that  G  (z)  is  non-increasing  in 
z.  Thus,  since  the  distribution  of  X^  is  stochastically  increasing 
in  A1  (A5) ,  we  conclude  that  I^X^X^)  =  E^  G  (X  )  is  non¬ 
increasing  in  A^,  The  same  argument  applies  to  the  other  components 
of  A  establishing  the  desired  result.  Q 

Suppose  that  A  is  the  non-negative  orthant  of  and  let  g 

be  the  simplex  g  =  {'A  :  E^=1  Ai=l}.  If  A4  and  A5  are  satisfied,  then 
by  P5  we  observe  that  for  any  A  e  g  and  real  c,  PcA{x<.x>  is  continu¬ 
ous  and  non- increasing  in  c .  If  the  lower  bound  of  P  ,  (Xr£.x)  as 

CA  ' 

C  ■*  00  is  less  than  a  for  all  x  e  X  and  all  A  £  §,  then  for  each  x  e  X 
and  A  e  g  there  exists  a  smallest  number  b  =  b(x,A)  such  that 

Pb(x  A)A^X=^  =  a*  10:16  confidence  bound  t(x)  defined  by  D2  is  then 

given  by 

t(x)  =  sup  9(b(x,A)A)  .  (2.7) 

Aeg 

Now  b(x, A)  is  easily  computed  using  root-finding  techniques  so  that 
the  computation  of  t(x)  reduces  to  searching  over  g  for  the  maximum 
of  9(b(x,A)A).  Many  routines  are  available  for  implementing  such 
searches.  For  the  situation  described  in  Corollary  Cl,  the  value  of 
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b  such  that  P^{x=^,  x}  -  a  is  unique  and  (2.7)  is  a  computationally 
feasible  version  of  (1.5). 

All  of  the  above  results  apply  mutatis  mutandis  to  the  con¬ 
struction  of  lower  confidence  bounds  and  hence  confidence  intervals. 
Applications  to  the  reliability  of  coherent  systems  involving  the 
binomial  or  other  distributions  are  possible.  In  particular,  the 
above  discussion  applies  directly  to  the  binomial  case  under  the 

transformation  A±  =  -  logU-p^),  i  =  1,2 . k; 

ea>  -  zj=1  -  -  log  n*=1  u-pt). 

3.  Systems  With  k  =  2 

As  was  noted  in  Section  1,  if  the  system  has  effectively  only 
one  component  (e.g.,  when  all  sample  sizes  are  equal),  then  the  pro¬ 
blem  reduces  to  the  well-known  case  of  finding  an  upper  confidence 
bound  for  a  single  Poisson  parameter.  Then  in  the  notation  of 
Section  1,  0(A)  =  A  and  if  t(x)  is  the  confidence  bound  for  A,  the 
lower  confidence  bound  (1.6)  for  reliability  R  becomes 

r(x)  =  1  -  t(x)/n  ,  (3.1) 

where  n  is  the  (common)  sample  size. 

The  two  component  case  (k  =  2)  presents  all  of  the  difficul¬ 
ties  of  the  general  case.  The  principal  problem  is  to  generate  an 
ordering  of  the  sample  points  x  =  (x^x^  which  will  lead  to  a  "good" 
confidence  bound  t(x)  computed  using  (1.5)  or  (2.7).  Several  differ¬ 
ent  methods  have  been  considered  and  implemented  to  varying  extents 
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In  the  course  of  this  investigation.  These  methods  may  be  described 
briefly  as  follows: 

(i)  The  x's  are  ordered  according  to  the  values  of  the  func- 

/o  o 

tion  t(x)  =  a-jX^  +  a2x2  +  Za  alxl  +  a2x2  *  where  za  is  the 
upper  a-th  quantile  of  the  standard  normal  distribution, 

(ii)  The  x’s  are  ordered  according  to  the  values  of  the  approx¬ 
imate  confidence  bound  obtained  from  the  maximum  likeli¬ 
hood  ratio  statistic  (see  Section  4). 

(iii)  The  ordering  is  generated  sequentially  by  considering  at 
each  stage  the  group  of  points  which  are  not  yet  ordered 
but  could  be  adjoined  without  violating  the  natural  par¬ 
tial  ordering  (see  D3  of  Section  2) .  The  next  point  in 
the  ordering  is  then  chosen  to  be  the  "best"  member  of  the 
candidate  group,  i.e.,  the  point  producing  the  smallest 
value  of  t (x)  given  by  (1.5). 

(iv)  The  ordering  is  chosen  so  as  to  minimize  E  (t(X)}  for  some 

G 

suitable  prior  distribution  G  over  the  values  of  X. 

(v)  The  ordering  is  generated  sequentially  in  the  manner  of 

(iii)  above  except  that  at  each  stage  the  candidate  points 
for  the  next  two  steps  are  considered  as  pairs  and  the 
next  point  selected  is  the  first  step  component  which, 
together  with  the  best  available  point  for  the  second 
step,  produces  the  smallest  sum  for  the  two  values  of 
t(x).  Note  that  the  point  that  appears  to  be  "best"  two 

steps  ahead  may  not  actually  be  chosen  when  that  stage  is 
reached. 
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It  is  clear  that  none  of  these  methods  is  special  to  the 
case  k  =  2.  Method  (i),  based  on  the  function  t(x),  which  is  really 
a  maximum  likelihood  estimate  of  an  asymptotically  valid  confidence 
bound,  was  used  to  generate  tables  of  bounds  for  the  case  k  =  2  in 
Johns  (1975) .  Method  (ii)  does  not  improve  substantially  on 
Method  (i)  for  moderate  values  of  the  X^'s.  Methods  (i)  and  (ii)  are 
asymptotically  equivalent  when  at  least  one  X^  becomes  large  (see 
Johns  1975)  and  indeed  standard  maximum  likelihood  results  guarantee 
that  both  are  asymptotically  optimal.  Method  (iii)  discussed  in 
Johns  (1977)  was  found  to  be  a  substantial  improvement  on  (i)  in  the 
strong  sense  that  when  the  Method  (iii)  ordering  is  used  the  values 
of  t(x)  are  often  smaller  and  only  rarely  slightly  larger  than  the 
values  for  corresponding  x's  produced  by  Method  (i) . 

The  semi-Bayesian  approach  of  Method  (iv),  which  minimizes 
the  expected  length  of  the  confidence  interval,  is  the  only  one  of 
the  five  that  is  directly  motivated  by  optimality  considerations. 

The  bound  resulting  from  any  reasonable  prior  must  at  least  be 
admissible.  In  pursuing  this  approach  it  was  decided  in  the  spirit 
of  objectivity  and  in  the  hope  of  rapid  convergence  to  asymptotic 
optimality  to  choose  a  prior  distribution  leading  to  an  unconditional 
probability  mass  function  for  ^  and  x2  constant  for  constant  values 
of  the  maximum  likelihood  estimator  a.jX^  +  a2x2  for  In 

particular,  the  prior  density  for  and  X^  was  taken  to  be 

-b  X  -b  A 

g(W  =  blb2  e  >  VX2  >  0  »  (3.2) 
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This  produces  the 


-3a1  -Ba 

where  =  (1-e  ) ,  =  (1-e  ) ,  B  >  0. 

unconditional  probability  mass  function 

-B(a1x1  +a2x2) 

p (x-^ * x2^  —  ^2.^2  ^  5 x2  —  *  •  ■>  •  (3.3) 

In  the  limiting  case,  as  B  +  0,  p(x^,x^)  becomes  essentially  uniform 
over  any  finite  set  of  points  (x^jX^) . 

The  actual  minimization  of  E  (t(X)}  may,  in  principle,  be 

Or 

accomplished  by  finding  the  ordering  which  minimizes  the  contribution 
to  Eg  among  all  orderings  of  length  N  where  N  may  be  arbitrarily 
large.  This  may  be  done  systematically  by  starting  at  the  origin 
(0,0)  and  constructing  a  tree  whose  nodes  at  each  stage  are  charac¬ 
terized  by  a  candidate  point  newly  adjoined  to  the  ordering  and  the 
corresponding  value  of  Ep(x)t(x),  where  the  sum  is  taken  over  all  x's 
occurring  in  the  path  leading  to  the  node,  including  the  one  just 
adjoined.  At  the  N-th  stage  the  node  having  the  smallest  accumulated 
sum  identifies  the  optimal  ordering  of  length  N.  This  process  may  be 
facilitated  by  eliminating  duplicate  nodes  and  discontinuing  branches 
when  a  node  is  reached  whose  value  exceeds  that  known  to  be  attain¬ 
able  in  N  stages.  Nevertheless,  because  of  the  rapid  increase  in  the 
number  of  nodes  considered  per  stage,  only  the  first  forty  or  so 
points  in  the  optimal  orderings  could  be  determined  even  using  a  very 
large  computer  facility. 

In  order  to  obtain  examples  of  admissible  orderings  with 
which  to  compare  the  results  of  other  methods,  this  computation  was 
performed  for  two  cases  using  an  IBM  370/168.  For  both  cases  the 
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values  a  =  .10  and  a^  =  .30  were  used.  For  the  first  case  the  proba¬ 
bilities  given  in  (3.3)  with  3=1  were  used  and  the  first  41  points 
of  the  optimal  ordering  were  obtained.  The  forty-first  stage  of  the 
computation  produced  4557  nodes.  For  the  second  case  the  limiting 
situation  as  3  0  where  the  p(x)'s  are  all  equal  was  used  and  the 

first  43  points  of  the  optimal  ordering  were  obtained.  The  number  of 
nodes  produced  at  the  forty-third  stage  was  6478. 

A  comparison  of  these  results  with  the  corresponding  results 
obtained  using  Methods  (i),  (iii),  and  (v)  is  indicated  in  Figure  1. 
The  horizontal  axis  indexes  the  first  50  points  in  the  ordering  pro¬ 
duced  by  Method  (v) ,  the  two-step  prospective  sequential  procedure. 
The  values  of  t(x)  for  these  indexed  points  for  the  five  methods  are 
indicated  by  the  plotted  symbols.  Values  of  t(x)  for  methods  other 

I 

than  (v)  are  shown  only  when  they  differ  from  those  produced  by  that 
method.  Based  on  this  evidence  it  appears  that  Method  (iv)  and 
Method  (v)  differ  very  little  and  that  both  are  better  than  the  other 
methods.  In  fact.  Method  (iv)  for  the  uniform  case  (3  =  0)  differs 
only  trivially  from  Method  (v) .  Since  the  use  of  Method  (iv)  for  the 
construction  of  tables  is  now  and  probably  always  will  be  impractical, 
we  are  led  to  the  choice  for  this  purpose  of  the  more  tractible  and 
virtually  equivalent  Method  (v) .  Prospective  sequential  methods 
looking  ahead  more  than  two  steps  might  be  feasible,  although  the 
complexity  of  the  computations  increases  rapidly  with  the  number  of 
steps.  However,  such  procedures  would  be  expected  to  produce  only 
minute  improvement  over  the  two-step  method. 
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Figure  1 .  Comparison  of  confidence  bounds  for  five  ordering  methods: 


Table  1  gives  values  of  the  upper  confidence  bound  t(x)  com¬ 
puted  using  (2.7)  with  the  sample  orderings  generated  by  the  two-step 
prospective  sequential  Method  (v) 0  The  values  are  given  for  the 
first  100  points  in  each  ordering  for  a  =  .01,  .05,  .10,  for 
al  =  .05 (.05) .45.  It  is  assumed  that  the  components  of  the  system 

are  indexed  so  that  a±  <  &2  which  Implies  a±  <  .50.  For  convenience 

) 

the  values  of  (x^»x2)  are  listed  systematically  rather  than  in  the 
order  generated  by  the  two-step  procedure.  This  table  provides  a 
basis  for  computing  accurate  confidence  bounds  for  the  case  k  =  2 
using  only  simple  interpolation.  If  values  of  aL  greater  than  .45 
but  (necessarily)  less  than  .50  are  required,  the  bound  for  a±  =  .5 
(corresponding  to  n^  =  n2)  may  be  used  for  interpolation.  This  bound 
is  obtained  by  simply  multiplying  the  ordinary  upper  confidence  bound 
for  a  single  Poisson  parameter  based  on  x  =  x1  +  x2  failures  by  .50 
(see,  e.g.,  Pearson  and  Hartley  1958  for  tables).  The  use  of  Table  1 
is  illustrated  by  the  following  two  examples: 

Example  1 .  Suppose  that  the  two  components  of  a  series  system  are 
tested  independently  using  sample  sizes  n^  =  300  and  n2  *  100 
respectively  with  the  corresponding  observed  numbers  of  failures 
Xx  =  3  and  X2  =  4.  Then  c  =  (1/300  +  1/100)  =  4/300  and 
al  -  l/cn^  =  .25  =  1  -  a2#  If  we  wish  to  find  a  95  percent  confi¬ 
dence  interval,  we  take  a  =  .05,  and  from  Table  1  we  find  the  confi¬ 
dence  bound  t(x)  for  9(A)  to  be  7.333.  Hence  by  (1.6)  the  95  percent 
lower  confidence  bound  for  system  reliability  R  is 
1  -  (4/300) (7.333)  =  .902. 
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Table  1.  The  Confidence  Bound  t(x)  for  k  =  2  for  the  First  100  Points  Generated  by  the 
Two-Stage  Optimal  Ordering  Method  for  Each  a.  and  a. 
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Table  1.  (Continued) 
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Table  1.  (Continued) 
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9.089 


A  somewhat  more  complicated  case  making  use  of  the  possibil¬ 
ity  of  combining  component  data  when  some  sample  sizes  are  equal,  as 
described  in  Section  1,  is  illustrated  by 

Example  2.  Consider  a  five-component  series  system  with  test  data 
=  0,  X2  =  4,  Xg  =  2,  X4  =  1,  X^  =  3,  based  on  corresponding  sam¬ 
ple  sizes  n^  =  n2  =  n3  =  300,  n^  =  n5  =  200.  Here 
c  =  (3/300  +  2/200)  =  .02,  =  a2  =  a3  =  l/300(.02)  =  1.6, 

=  a^  =  l/200(.02)  =  1/4.  For  the  equivalent  k*  =  2  problem  we 

divide  the  two  distinct  values  of  the  a. *s  by  their  sum  c  =  1/6  + 

i  o 

1/4  =  5/12,  obtaining  a*  =  .40,  a*  =  .60.  The  corresponding  numbers 
of  failures  are  X^  =  X^^  +  X2  +  X3  =  6  and  X*  =  X^  +  X5  =  4.  For  the 
reduced  problem  we  consult  Table  1  for  a  =  .10  to  find  the  90  per¬ 
cent  confidence  bound  7,564.  This  must  be  multiplied  by  cq  to 
obtain  the  bound  for  the  original  0(X).  The  90  percent  lower  bound 
r(x)  for  the  system  reliability  given  by  (1.6)  thus  becomes 
1  -  ccQt(x)  =  1  -  (.02) (5/12) (7.564)  =  .937. 

If  the  observed  numbers  of  failures  fall  outside  the  limits 
of  Table  1,  the  approximate  methods  discussed  in  the  next  section 
may  be  employed. 
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4 .  Approximations 

In  this  section  approximate  confidence  bounds  for  systems 
with  more  than  two  components  having  distinct  sample  sizes  are 
developed.  The  use  of  the  maximum  likelihood  ratio  confidence  bound 
for  cases  falling  outside  the  scope  of  Table  1  is  also  discussed. 

For  k  >  2  it  is  impractical  to  generate  the  optimal  ordering 
and  the  corresponding  values  of  the  upper  confidence  bound  for  more 
than  a  few  illustrative  cases.  Thus,  some  method  of  approximating 
solutions  for  k  >  2  with  acceptable  precision  is  required.  The 
approach  to  be  followed  here  is  to  find  a  k  =  2  problem  which  is  suf¬ 
ficiently  similar  in  structure  to  the  given  k  >  2  problem  so  that  the 
confidence  bounds  for  the  two  problems  are  essentially  the  same  except 
for  a  normalizing  factor.  This  method  may  be  thought  of  as  an  exten¬ 
sion  and  refinement  of  the  Lindstrom-Madden  procedure  (see  Lloyd  and 
Lipow  1977,  and  cf .  Harris  and  Soms  1980). 

The  Lindstrom-Madden  method  first  estimates  the  reliability 
by  maximum  likelihood  and  then  uses  the  k  =*  1  confidence  bound  solu¬ 
tion  for  the  component  with  the  smallest  sample  size  and  a  fictitious 
number  of  failures  determined  so  as  to  reproduce  the  estimated  system 
reliability.  The  procedure  proposed  here  is  to  estimate  two 
quantities,  the  value  of  0  and  the  variance  of  the  maximum  likelihood 
estimate  of  0,  and  to  use  these  estimates  to  find  a  k  *»  2  problem 
based  on  two  of  the  original  a^'s  with  a  pair  of  corresponding  ficti¬ 
tious  observation  values  chosen  to  reproduce  the  estimated  quantities . 
The  two  a^’s  are  chosen  to  be  as  large  as  possible  (corresponding  to 
sample  sizes  as  small  as  possible)  subject  to  two  constraints.  The 
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first  constraint  guarantees  that- the  resulting  fictitious  observations 
are  non-negative.  The  second  constraint  requires  that  a^’s  corre¬ 
sponding  to  zero  failures  in  the  original  problem  not  be  considered 
unless  all  but  one  of  the  components  have  zero  failures .  These  con¬ 
siderations  lead  to  a  unique  k  =  2  problem  whose  solution  provides  a 
very  good  approximation  to  the  solution  for  the  original  k>  2  problem. 

The  elimination  of  a^’s  for  components  exhibiting  zero 
failures  is  justified  by  the  fact  that  the  maximum  likelihood  ratio 
confidence  bound  (discussed  later)  is  invariant  under  such  transfor¬ 
mations.  That  is,  if  the  dimension  k  is  reduced  by  the  elimination 
of  all  a^'s  corresponding  to  zero  failures,  then  the  value  of  the 
maximum  likelihood  ratio  bound  remains  unchanged. 

The  choice  of  the  two  quantities  whose  estimates  are  used  to 

determine  the  pair  of  fictitious  observations  is  supported  by  analogy 

with  Lindstrom-Madden  (in  the  case  of  0)  and  by  the  fact  that  the  two 

estimates  are  the  ingredients  of  the  asymptotic  maximum  likelihood 

ratio  confidence  bound  thereby  insuring  asymptotic  optimality.  The 

details  of  the  approximation  algorithm  are  as  follows: 

ic 

(a)  First  0(A)  =2.  .  a. A.  is  estimated  by 

i=l  ii 

k 

0=2  a.X.  ,  (4.1) 

i=l  1  1 

a  k  2 

and  the  quantity  Var(0)  =  2^_^  a^A^  *-s  estimated  by 

a  k  ? 

v  =  2  al  X.  .  (4.2) 

i  i 
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(b)  Next  if  at  least  one  component  exhibits  one  or  more  failures 


the  pair  (a^a.),  i  <  j,  is  selected  so  that 

/s 

(i)  a^  _<  v/0  ,  and 


(ii)  a^  >_  v/0 

Subject  to  these  conditions,  a^  and  a^.  are  taken  to  be  the 

largest  available  values  associated  with  at  least  one 

failure.  If  all  a^'s  satisfying  (i)  correspond  to  zero 

failures,  then  a^  is  taken  to  be  the  largest  in  that  group. 

If  all  the  a_.  's  satisfying  (ii)  correspond  to  zero  failures, 

then  a.  is  taken  to  be  a^.  If  all  components  exhibit  zero 

failures,  v/0  is  indeterminate  and  a.  and  a.  are  taken  to  be 

1  3 

a^  and  a^  respectively.1 

(c)  The  pseudo-observations  x*  and  x*  are  computed  by  the 
formulae 


a .0 -  v 
JL 


x*  = 

1  a.  (a.  -  a  ) 
1  3  i 


(4.3) 


= 


v-  a^0 


These  values  will  be  non-negative  by  conditions  (i)  and  (ii) 


of  (b) ,  and  when  associated  with  a^ 


and  a. 
3 


respectively  they 


1For  this  case  the  resulting  confidence  bound  is  exact  and  is  the 
same  as  would  be  obtained  by  multiplying  a^  times  the  upper  confi¬ 
dence  bound  for  a  single  Poisson  parameter  when  zero  failures  are 
observed,  i.e.,  t(0)  =  a^  log  (1/a) . 


reproduce  the  values  of  9  and  v  provided  all  other  observa¬ 
tions  are  replaced  by  zeros. 

(d)  The  k  =  2  problem  with  and  x*  associated  with 

a*  =  a.j/ (a^+  aj)  and  a*  =  a^/Ca^a^)  respectively  may  now 
be  treated  using  Table  1  to  yield  t(x*,x*).  Since  x*  and  x* 
are  not  necessarily  integers,  it  may  be  necessary  to  inter¬ 
polate  with  respect  to  these  arguments  as  well  as  the  value 
of  a^  =  a* . 

(e)  The  approximate  upper  confidence  bound  t*  for  9  for  the 
original  k  >  2  problem  is  then  given  by 


t*(x)  =  (a±+ aj)t(x*,x*) 


(4.4) 


In  order  to  check  the  validity  of  this  approximation 
algorithm,  the  confidence  bounds  for  the  first  24  points  in  the  opti¬ 
mal  (two-stage  prospective)  ordering  were  computed  for  a  typical  k=  3 
case.  These  results  were  obtained  by  a  rather  laborious  method  based 
on  formula  (2.7)  and  involving  repeated  interactive  searches  of  the 
(X^,X2»'>13)  simplex.  The  approximation  algorithm  was  applied  to  each 
of  these  points  and  the  comparative  results  are  shown  in  Table  2. 

The  values  of  a^,  a.^,  and  a^  for  this  example  were  chosen  to  be  of 
roughly  the  same  magnitude  but  not  so  close  that  the  combination  of 
any  two  would  be  indicated. 

To  check  the  algorithm  for  cases  farther  from  the  origin 
several  additional  examples  were  considered  using  a  somewhat  differ¬ 
ent  method  which  avoids  the  necessity  for  sequentially  generating  the 
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Table  2.  The  Performance  of  the  Approximation  Algorithm  for  the 
First  24  Ordered  Points  for  a  Typical  Example  With  k  =  3. 


a  =  .05;  a^  =  .2,  a^  =  .3,  a^  =  .5 


Two-Step  Prospective 
Sequential 


Approximation 

Algorithm 


Relative 


n 

X1 

X2 

X3 

t(x) 

t*(x) 

Error 

1 

0 

0 

0 

1.498 

1.498 

+  00.0  % 

2 

1 

0 

0 

1.553 

1.555 

+  00.1 

3 

0 

1 

0 

1.656 

1.663 

+  00.4 

4 

2 

0 

0 

1.703 

1.705 

+  00.2 

5 

1 

1 

0 

1.756 

1.540 

-  12.3 

6 

2 

1 

0 

1.861 

1.752 

-  05.9 

7 

3 

0 

0 

1.933 

1.891 

-  02.1 

8 

0 

2 

0 

1.981 

1.995 

+  00.7 

9 

1 

2 

0 

2.085 

2.052 

-  01.6 

10 

4 

0 

0 

2.088 

2.094 

+  00.3 

11 

3 

1 

0 

2.162 

1.991 

-  07.9 

12 

2 

2 

0 

2.268 

2.208 

-  02.7 

13 

5 

0 

0 

2.300 

2.309 

+  00.4 

14 

0 

3 

0 

2.354 

2.397 

+  01.8 

15 

0 

0 

1 

2.372 

2.372 

+  00.0 

16 

4 

1 

0 

2.382 

2.291 

-  03.9 

17 

1 

0 

1 

2.429 

2.431 

+  00.1 

18 

1 

3 

0 

2.472 

2.438 

-  01.4 

19 

5 

1 

0 

2.504 

2.486 

-  00.7 

20 

0 

1 

1 

2.543 

2.529 

-  00.5 

21 

6 

0 

0 

2.575 

2.551 

-  01.0 

22 

2 

0 

1 

2.602 

2.589 

-  00.8 

23 

2 

3 

0 

2.609 

2.641 

+  01.2 

24 

1 

1 

1 

2.660 

2.668 

+  00.3 
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optimal  ordering  of  the  sample  points.  For  these  examples  the  order¬ 
ing  was  generated  by  the  values  of  the  maximum  likelihood  ratio  con¬ 
fidence  bounds  (see  below)  associated  with  the  sample  points.  Since 
this  ordering  is  asymptotically  optimal  (for  large  X^’s),  it  can  be 
expected  to  produce  good  results  for  sample  points  well  removed  from 
the  origin.  The  values  of  the  confidence  bounds  calculated  using 
(2.7)  and  the  approximations  obtained  by  the  proposed  algorithm  are 
shown  for  these  examples  in  Table  3. 

In  Table  3  the  tendency  of  the  computed  values  of  t(x)  to 
be  slightly  larger  than  the  algorithm  values  may  be  due  to  the  fact 
that  the  ordering  used  to  compute  the  former  is  non-optimal.  Formal 
application  of  the  algorithm  may  occasionally  result  in  the  selection 
of  nearly  equal  values  a^  and  a_.  in  step  (b)  .  When  this  happens, 
improved  results  may  be  obtained  by  first  reducing  the  dimension  k  by 
combining  the  nearly  equal  a^'s  and  then  applying  the  algorithm. 

This  method  was  used  for  the  two  cases  in  Table  3  marked  by  (+)  . 

The  algorithm  values  for  all  cases  where  a  =  ( .14, .16, .70)  or 
a  =  ( .15, .41,  .44)  do  not  differ  substantially  from  the  values  which 
would  be  obtained  by  reducing  to  the  k  =*  2  case  by  combining  the 
nearly  equal  a^'s.  Similarly,  the  values  for  the  cases  where 
a  =  ( .32, .33, .35)  can  be  nearly  reproduced  by  multiplying  the  k  =  1 
bound  by  a  3  .333. 

The  number  of  sample  points  appearing  in  the  ordering  before 
a  given  level  of  t(x)  is  reached  increases  rapidly  as  k  increases. 
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Table  3.  Several  Examples  Comparing  the  Algorithm  Values  With  the 
Exact  Bounds  Based  on  the  Ordering  Generated  by  the  m£rb. 


a  =  .05 


al  a2  a3 

X1  X2 

X3 

Position  in 
Ordering 

t(x) 

Algorithm 
Value  t*(x) 

m£rb 

.20, .30,  .50 

2,  2, 

1 

50 

3.329 

3.273 

3.070 

5,  0, 

2 

100 

4.068 

3.957 

3.798 

6,  6, 

0 

200 

4.887 

4.789 

4.708 

9,  1, 

3 

400 

5.980 

5.892 

5.795 

.14, .16,  .70 

5,  4, 

0 

50 

3.026 

3.084+ 

2.218 

2,10, 

0 

100 

3.424 

3.531+ 

2.921 

2,  5, 

1 

200 

4.407 

4.086 

3.729 

.15,. 41, .44 

13,  0, 

0 

50 

3.245 

3.223 

2.980 

5,  2, 

1 

100 

4.028 

3.860 

3.699 

3,  3, 

2 

200 

4.861 

4.738 

4.559 

.32, .33, .35 

0,  3, 

2 

50 

3.492 

3.518 

3.256 

i,  5, 

1 

100 

4.339 

4.332 

4.076 

1,  5, 

3 

200 

5.198 

5.234 

4.993 

Reduction  to  k=2  case  by  combining  a ^  and  a^. 


38 


Hence  for  k  >  2  the  algorithm  may  be  applied  to  sample  points  posi¬ 
tioned  quite  far  along  in  the  ordering  without  requiring  values 
beyond  the  scope  of  Table  1,  as  is  seen  in  the  examples  of  Table  3. 

The  results  detailed  in  Tables  2  and  3  suggest  that  the 
application  of  the  proposed  algorithm  together  with  the  combining  of 
nearly  equal  a^'s  when  indicated  will  nearly  always  produce  confi¬ 
dence  bounds  for  0 (X)  subject  to  relative  errors  not  exceeding  10 
percent  and  often  much  less.  Furthermore,  the  lower  confidence  bound 
for  reliability  r(x)  given  by  (1.6)  will  exhibit  a  much  smaller  rela¬ 
tive  error  since  the  relative  error  in  approximations  for  t(x) 
applies  only  to  the  difference  between  the  lower  bound  and  one;  a 
quantity  which  at  worst  is  of  the  order  of  1/10  in  the  contemplated 
applications . 

In  situations  where  very  large  sample  sizes  are  available,  it 
may  happen  that  the  observed  numbers  of  failures  exceed  the  limits  of 
Table  1  even  though  the  system  reliability  is  high.  For  such  cases 
the  maximum  likelihood  ratio  bound  (m£rb)  may  be  used.  This  approxi¬ 
mate  confidence  bound  is  obtained  in  the  usual  way  from  the  maximum 
likelihood  ratio  statistic  for  testing  the  hypothesis  :  0(A)  =  0^ 
versus  all  alternatives.  The  corresponding  one-sided  confidence 
bound  may  be  shown  (see  Johns  1975)  to  be  determined  as  follows:  Let 
y  be  the  positive  real  root  less  than  1/ (largest  a^  for  which  >  0) 
of  the  equation 
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(4.5) 


XU  ’  2  £  Xi 


a . A 
_ ILL 


1-  a^ii 


+  log  (1  —  a_jii) 


2 

where  X-^  2a  :*’s  t*ie  uPPer  100-(2a)-th  percentage  point  of  the  chi- 
squared  distribution  with  one  degree  of  freedom.  Then  the  quantity 

k 

t(X)  =  Z  a.X  / (1-  a.y)  (4.6) 

i=l  1  1  1 

is  the  approximate  upper  1-  a  level  confidence  bound  for  ©(A^)  .  As 

A 

max(A- . .  .  ,A  )  °°,  the  m&rb  t(X)  may  be  shown  (see  Johns  1975)  to 

12  n 

be  asymptotically  equivalent  to 


t(X) 


k 

Z 

i=l 


a.X.  +  z 
i  i  a 


k 

Z 

i=l 


a2X. 


(4.7) 


where  z^  is  the  lOOa-th  percentage  point  of  the  standard  normal 

distribution. 

Neither  of  these  approximate  bounds  is  useful  for  sample 
points  near  the  origin  in  the  usual  orderings,  but  t  given  by  (4.6) 
becomes  sufficiently  precise  for  application  to  sample  points  beyond 
the  scope  of  Table  1.  A  comparison  of  t(x^,X2)  and  the  corresponding 
m&rb  for  the  last  (i.e.,  the  100-th)  points  in  each  of  the  optimal 
orderings  for  the  cases  covered  in  Table  1  is  given  in  Table  4. 
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Table  4.  m£rb/t(x1,x2)  for  the  100-th  Point  in  Each  Ordering 


al 


a 

.05 

.10 

.15 

.20 

.25 

.30 

.35 

.40 

.45 

.01 

.62 

.95 

.97 

.97 

.98 

.97 

.98 

.97 

.99 

.05 

.72 

.88 

.96 

.95 

.98 

.97 

.99 

.95 

.96 

.10 

.89 

.95 

.95 

.96 

.97 

1.00 

.97 

.96 

.98 

These  results  suggest  that  the  m&rb  possesses  satisfactory 
precision  for  sample  points  beyond  those  listed  in  Table  1  for  k  =  2 
whenever  a^  _>  .10.  The  last  column  of  Table  3  giving  the  m&rb  values 
for  the  examples  considered  illustrates  the  fact  that  for  k  >  2  the 
m£rb  tends  to  underestimate  the  correct  value  of  the  bound  for  sample 
points  within  the  range  that  can  be  dealt  with  using  Table  1  and  the 
algorithm . 

Two  potential  sources  of  error  for  the  lower  confidence  bound  on 

system  reliability  remain  to  be  discussed.  They  are  (i)  the  Poisson 

approximation  to  the  binomial  distribution  of  the  observed  failures,  and 

(ii)  the  approximation  for  reliability  given  in  (1.2)  and  reflected  in 

the  formula  (1.6)  for  the  lower  bound  r(x) .  It  is  intuitively  clear  from 

(1.1)  f f .  that  the  "worst  case"  for  the  Poisson  approximation  should 

occur  when  k  =  1;  since  to  match  a  given  k  =  1  level  of  reliability,  say 

1-p,  by  a  k  >  1  case,  we  must  have  p  =  .  p . ,  so  that  the  p.'s  must  be 

i=l  i  r i 

smaller  than  p  which  tends  to  improve  the  Poisson  approximation. 
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For  the  case  k  =  1  the  familiar  upper  confidence  bounds  for  a 
single  Poisson  parameter  apply  and  the  actual  coverage  probabilities  for 
the  proposed  method  (3.1)  can  be  computed  for  any  n  and  p  from  tables  of 
the  binomial  distribution.  The  results  of  several  such  calculations  are 
shown  in  Table  5. 


Table  5.  "Worst 
for  r(x) . 

Case"  Analysis  (k=  1)  . 

Minimum  Coverage 

Probabilities 

Reliability 

1  -  a 

q  -  i  -  P 

.90 

.95 

.99 

.95 

.906 

.954 

.991 

.90 

.912 

.958 

.992 

.80 

.924 

.965 

.995 

.70 

.936 

.972 

.997 

These  values  suggest  that  the  approximations  operate  to  make  the 
proposed  confidence  bounds  slightly  conservative.  It  is  interesting  to 
observe  that  the  minimum  coverage  probabilities  are  not  drastically  dif¬ 
ferent  from  the  nominal  values,  even  for  a  true  reliability  as  low  as 
.70. 
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